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LONG-TERM  GOALS 

The  goal  of  this  project  is  to  improve  our  understanding  of  turbulence  and  small  scale  processes  in  the 
aquatic  near  surface  layer  and  their  relation  to  surface  waves  and  meteorological  forcing.  An  improved 
understanding  of  these  processes  should  result  in  the  improvement  of  turbulence  parameterization 
schemes  used  in  the  aquatic  surface  layer  and,  as  a  result,  in  more  accurate  model  predictions. 

OBJECTIVES 

Our  approach  is  to  synthesize  the  available  near-surface  turbulence  kinetic  energy  (TKE)  dissipation 
data  sets,  along  with  their  simultaneous  meteorological,  surface  wave,  and  current  measurements 
(when  available).  Then  classify  how  various  surface  forcing  conditions  affect  the  results.  Specific 
objectives  of  this  analysis  are  to  provide: 

1)  A  consistent  assessment  of  existing  data  sets. 

2)  Recipe/s  of  parameterization  schemes  of  TKE  dissipation  rates,  s,  as  a  function  of  the  forcing  (i.e. 
wind  stress,  wave  age,  wave  height,  buoyancy  flux,  stratification).  Comparison  to  one-dimensional 
model  runs,  including  surface  wave  effects,  will  help  to  reach  this  objective  and  may  later  be 
incorporated  into  three  dimensional  models. 

3)  Estimates  of  the  fraction  of  the  surface  energy  flux  into  the  ocean  (via  the  surface  waves)  which  is 
parted  to  the  mean  and  the  turbulent  flow. 

4)  Guidance  for  the  design  of  future  experiments  to  fill  in  currently  existing  gaps  in  our  knowledge. 

APPROACH 

We  have  collected  data  sets  from  various  investigators,  who  agreed  to  contribute  relevant  data,  as  well 
as  from  the  published  literature.  Data  includes  near  surface  TKE  dissipation  rates  measured  from  the 
following  platforms:  free  rising/falling  turbulence  profilers,  quasi-horizontal  gliders  equipped  with 
shear  probes  and  fast  thermistors,  submarine  and  ship-bow  mounted  turbulence  sensors,  acoustic 
travel-time  current  meters,  and  drag  spheres.  In  some  cases  hydrographic  data,  comprising  vertical 
profiles  of  temperature  and  salinity,  and  hence  density,  were  available.  Atmospheric  boundary  layer 
(ABL)  fluxes  (wind  stress  and  heat  flux)  were  available  in  some  instances  while  for  most  only  wind 
speed  or  surface  stress  data  were  available.  For  most  data  sets  various  measured  wave  parameters  such 
as  significant  wave  height,  peak  frequency  (or  period,  or  phase  speed)  were  available  while  for  the 
others  only  wave  heights  were  available.  Data  were  collected  both  in  surface  layers  (SL)  of  the  ocean 
and  lakes,  at  various  geographical  locations,  and  under  a  variety  of  atmospheric  and  surface  wave 
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conditions.  The  latter  is  of  importance  since  a  variety  of  studies  have  shown  a  correlation  between 
wave  age  and  the  rate  of  energy  input  from  the  wind  to  waves  and  hence  to  turbulence  (e.g.  Terray  et 
ah,  1997). 

The  available  TKE  dissipation  rates,  hydrographic,  meteorological,  wave,  and  data  have  been  put  into 
a  database  and  Matlab  programs  were  written  to  facilitate  analyses  and  visualization  of  the  data  (see 
summary  of  available  data  in  Table  1).  The  datasets  facilitate  statistical  analysis  of  various  parameters 
such  as  the  distribution  of  TKE  dissipation  rates  as  a  function  of  wind  stress,  wave  height  and  age, 
stratification  and  stability  in  the  oceanic  boundary  layer  (OBL),  etc. 

Identification  of  the  dominant  forcing  parameters  (e.g.  surface  friction  velocity,  wave  age, 
absence/presence  of  wind  waves,  absence/presence  of  swell)  for  each  of  the  dissipation  data  sets 
allows  the  detennination  of  the  best  apparent  nondimensional  parameterization  groups  on  which  to 
scale  the  dissipation  in  the  wind-wave  zone  and  the  layer  beneath  it.  Estimates  of  the  dominant 
turbulence  length  scales  and  time  scales  involved  in  the  near  surface  turbulence  processes  can  be  made 
and  related  to  the  observed  dissipation  rates  and  the  pertinent  forcing.  For  cases  of  enhanced  TKE 
dissipation  the  datasets  can  be  used  to  determine  at  what  depth  range,  the  commonly  used  constant 
stress  layer  scaling  and/or  convective  scaling  become  valid. 

Table  1.  Summary  of  available  meteorological  and  wave  data  for  the  various  experiments  analyzed 
(ranges  are  given  in  most  cases):  U10is  the  10-  m  wind  speed,  z is  the  surface  wind  stress,  u  wis  the 
water  surface  friction  velocity,  Hs  is  the  significant  wave  height,  Tp  is  the  period  at  the  peak  of  the 
wave  spectrum,  and  F  is  the  rate  of  wind  energy  input  into  the  waves. 


Instrument 

U,0 

[m/s] 

T 

[N/m2] 

tl 

[m/s] 

Hs 

[m] 

TP 

[m] 

Fx  104 
[m3/s3] 

Refer. 

Hot  film 

anemometer 

6.0 

0.051 

0.007 

0.40 

1.80 

0.65 

Stewart  and 
Grant  (1962) 

Miniature 

dragsphere 

5.8-11.2 

0.069-0.237 

0.008-0.015 

0.26-58 

1. 7-3.3 

0.65-3.23 

Kitagorodskii 
et  al.  (1983) 

Shear  probes  (free 
falling  profiler) 

8.6 

0.104 

0.010 

3.3  (swell) 

15.0  (swell) 

6.94 

Gregg  (1987) 

Shear  probes 
(submarine) 

5 .2-8.2 

0.041-0.102 

0.006-0.010 

0.56 

4.0 

0.77-2.97 

Osborn 

(1992) 

Shear  probes  (free 
rising  profiler) 

12.6 

0.252 

0.016 

1.0  (wind) 

3.0  (swell) 

4.0  (wind) 

12.0  (swell) 

12.7 

Anis  and 
Mourn  (1995) 

BASS  and 
Dragsphere 

6.9-15.7 

0.104-0.902 

0.010-0.030 

0.16-0.49 

1. 4-2.4 

0.96-18.2 

Terray  et  al. 
(1996) 

Acoustic  current 
meter  (bow 
mounted) 

8.0-11.8 

0.079-0.377 

0.009-0.019 

0.88-2.62 

3.4-6. 3 

0.78-6.42 

Drennan  et  al. 
(1996) 

Shear  probes  (quasi¬ 
horizontal  glider) 

6. 9-9. 9 

0.056-0.127 

0.005-0.012 

0.62-1.25 

3. 6-4.3 

1.53-14.2 

Greenan  et  al. 
(2001) 

Shear  probes  (free 
falling  profiler) 

6.8-10.2 

0.046-0.194 

0.005-0.012 

0.55-1.61 

3. 6-4.3 

2.20-10.3 

Greenan  et  al. 
(2001) 

EM  velocity  probes 
(bow) 

16.8-19.2 

0.067-0.779 

0.008-0.028 

0.97-4.32 

4.5-11.6 

2.30-56.5 

Soloviev  and 
Lukas  (2003) 

Based  on  the  wind  and  wave  data,  the  fraction  of  the  wind  energy  flux  Ejq  in  the  ABL  that  is  dissipated 
in  the  OBL  under  various  wind/wave-age/wave-height  conditions  is  quantified  (These  estimates  range 
between,  roughly,  1%  to  10%;  e.g.  Anis  and  Mourn,  1995).  Currently  researchers  are  in  disagreement 
as  to  the  amount  of  energy  flux  from  the  waves  to  the  ocean  (i.e.  into  the  surface  currents).  For 
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example  Crawford  and  Large  (1996)  assume  that  only  a  negligible  amount  of  the  energy  that  enters  the 
wave  layer  indeed  goes  into  the  ocean  currents.  Although  this  assumption  may  work  well  for  climate 
models  it  disagrees  with  results  from  several  experimental  and  theoretical  studies  and  may  be  an 
inaccurate  assumption  for  forecast  models  of  currents  and  waves.  More  importantly,  if  indeed  10%  of 
Eio  is  dissipated  in  the  OBL  this  may  well  be  a  manifestation  that  a  non-negligible  amount  of  energy 
flux  enters  the  ocean  and  goes  either  into  the  mean  current  field  or  into  turbulence  and  in  both  cases 
TKE  dissipation  rates  may  be  enhanced. 

WORK  COMPLETED 

Most  of  the  available  data  sets  from  various  field  experiments  during  which  TKE  dissipation  rates  and 
wave  measurements  were  carried  out  in  the  near  SL  of  ocean  and  lakes  have  been  acquired  (see  Table 
1  for  a  summary).  A  variety  of  instruments  and  platforms  have  been  used  to  measure  the  TKE 
dissipation  rates  and  while  several  of  the  data  sets  include  detailed  surface  wave  measurements,  some 
had  only  wave  estimates  (height  and  period)  from  the  ship's  bridge  observations.  Wind  speeds  for  the 
various  data  sets  that  were  analyzed  ranged  from  5.2-19.2  m/s,  and  wave  heights  and  periods  were 
between  0.16-4.32  m  and  1.4-15.0  s,  respectively  (see  Table  1).  Several  of  the  TKE  dissipation  data 
sets  include  relatively  shallow  measurements  on  the  order  of  a  few  significant  wave  heights  (Drennan 
et  ah,  1996;  Terray  et  ah,  1996;  Soloviev,  2003),  while  some  of  the  data  sets  include  deeper 
measurements  (e.g.  Anis  and  Mourn,  1995;  Greenan  et  al  2001;  Gregg  1987)  spanning  the  OBL. 

First,  a  comparison  of  the  various  data  sets  was  carried  out  using  a  common  framework  of  scaling  TKE 
dissipation  estimates  in  wall-layer  coordinates  in  which  the  dimensionless  TKE  dissipation  rate  is 
given  by  s/(u*w  /kz)  and  the  dimensionless  depth  by  gz/u*w  .  If  wall-layer  scaling  is  indeed  valid  we 
expect  that  s/(u*3/kz )  ~  1.  Fig.  1  shows  this  scaling  for  6  of  the  datasets  of  Table  1. 


Figure  1.  Dimensionless  dissipation  rate,  s/(u3/kz),  as  a  function  of  dimensionless  depth,  gz/u*2. 

The  constant  stress  layer  is  represented  by  the  vertical  line  (red),  s/(u*3/kz)  =  1  (sis  the  TKE 
dissipation  rate,  u  is  the  friction  velocity  in  water,  k=  0.4  is  von  Karman  's  constant,  and  z  is  the 
depth).  The  upper  panels  represent  dissipation  estimates  taken  in  relatively  shallow  depths  (a  few 
meters  at  most)  from  ship-bow  mounted  sensors  and  from  moored  instruments.  The  lower  panels 
represent  dissipation  estimates  made  from  profiling  instruments,  which  on  average  extend  to  greater 
depths  (10  m  or  more).  Sources  for  the  data  sets  are  noted  on  the  individual  panels. 
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Next,  we  examined  the  statistical  distribution  of  the  scaled  dissipations,  which  showed  that  the  wall 
layer  scaling  severely  underestimated  TKE  dissipation  in  the  SL.  More  specifically,  bootstrap  statistics 
of  the  wall  layer  scaling  show  that  mean  values  of  s/(u*3/kz)  range  from  3.6  and  up  to  34.5  (Table  2). 

A  different  test  of  the  validity  of  wall-layer  scaling  was  carried  out  by  examining  the  fraction  of  wind 
energy  flux  in  the  ABL,  Eio,  which  is  dissipated  in  the  OBL  (Eio=  xUio ;  x  is  the  surface  wind  stress 
and  Uio  is  the  wind  speed  at  10  m  height).  Wall  layer  scaling  predicts  the  dissipation  to  be  ~1%  of  Eio 
(Oakey  and  Elliott,  1982).  Laboratory  and  field  measurements  in  combination  with  a  model  used  by 
Richman  and  Garrett  (1977)  predicted  a  higher  percentage  (4-9%).  We  have  vertically  integrated 


Instrument 

Eio 

[W/m2] 

s/Eio 

r%i 

u*/kz  xlO7 

s/(u*/kz) 

Reference 

Shear  probes  (free 
falling  profiler) 

0.89 

14.6 

0.90-34.2 

30.2  (18.2-41.2) 

Gregg  (1987) 

Shear  probes 

(submarine) 

0.21-0.83 

5.97  (3.25-8.03) 

0.63-49.4 

3.6  (2. 6-4. 8) 

Osborn  (1992) 

Shear  probes  (free 
rising  profiler) 

3.16 

10.2 

7.2-64.5 

6.5  (3.1-10.8) 

Anis  and  Mourn  (1995) 

BASS  and  Dragsphere 

0.72-14.2 

1.94  (1.36-2.67) 

12.0-1230 

4.9  (3 .4-6.7) 

Terray  et  al.  (1996) 

Acoustic  current  meter 
(bow) 

0.66-3.89 

11.2  (7.92-15.1) 

9.79-104 

34.5  (23.3-50.0) 

Drennan  et  al.  (1996) 

Shear  probes  (free 
falling  profiler) 

0.18-1.29 

6.0  (3 .9-8.4) 

0.36-10.1 

9.7  (6.7-12.9) 

Greenan  et  al.  (2001) 

Shear  probes  (quasi¬ 
horizontal  glider) 

0.15-1.43 

4.7  (2. 5-7. 5) 

0.34-22.1 

11.0  (9.1-13.4) 

Greenan  et  al.  (2001) 

EM  velocity  probes 
(bow) 

0.45-14.9 

1.87  (1.84-1.89) 

4.96-657 

6.2  (6. 1-6.3) 

Soloviev  and  Lukas  (2003) 

Table  2.  Summary  of  scaling  results  for  various  experiments.  Variables  are  defined  in  the  text 
(bootstrap  mean  values  and  95%  confidence  intervals  are  given  in  parentheses). 

dissipation  rates,  si,  in  the  SL  (depths  on  the  order  of  a  few  wave  heights  and  less  than  10  m  for  all  sets 
in  this  analysis).  The  results  are  summarized  in  Table  2  and  show  that  in  the  presence  of  waves  the 
values  are  significantly  higher  than  the  wall  layer  predicted  1%  (one  has  to  consider  si  computed  here 
as  a  lower  limit  on  the  total  wind  energy  flux  dissipated  in  the  OBL  since  we  included  only  the  SL). 
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Figure  2.  Example  of  dissipation  rates  vs.  depth  in  the  scaled  wave  coordinates  of  Terray  et  al. 
(1996).  The  solid  line  (red)  represents  the  best  fit  of  their  WAVES  data:  eH/F  =  0.3(~/Hf2. 


Terray  et  al.  (1996)  suggested  that  in  the  presence  of  waves,  when  wall-layer  scaling  underestimates 
dissipation  in  the  upper  SL,  TKE  dissipation  will  be  dependent  on,  F,  the  rate  of  energy  input  from  the 
wind  to  the  waves  and,  Hs.  The  following  was  found  to  hold:  sHJF  =  0.3 (z/Hf  for  z  >0.6HS  and  sHJF 
=  0.83  for  z  <0.6HS.  Below  the  wave  affected  layer  the  “conventional”  wall-layer  may  be  valid.  We 
have  binned  the  scaled  dissipation  values,  sHJF  in  scaled  depth  bins,  z/ Hs,  and  computed  the  mean 
and  95%  bootstrap  confidence  intervals  for  each  bin,  using  the  values  of  F  and  Hs  given  in  Table  1 . 
Results  of  this  scaling  are  presented  in  Fig.  2. 


Profiles  of  TKE  dissipation  rates,  s,  were  modeled  using  the  1-  dimensional  turbulence  model  GOTM. 
The  scheme  used  was  the  so-called  generic  two-equation  model  (for  details  see  Umlauf  and  Burchard, 
2003).  Physically,  the  flow  beneath  the  surface  in  presence  of  waves  is  taken  as  a  generalization  of 
pure  shear-free  turbulence.  In  our  model  runs  we  have  “injected”  TKE  from  the  waves  at  several 
depths  levels,  z0,  expressed  in  terms  of  the  significant  wave  height,  e.g.  z0=  0.25//  to  zq=  1.0//,.  The 
injection  level  might  be  regarded  as  the  so-called  “breaking  depth”  (Terray  et  al.,  1996).  The 
magnitude  of  the  TKE  flux  injected  was  taken  simply  as  a  constant  factor,  a,  times  u*w  ,  and  spatial 
(vertical)  TKE  decay  rates  chosen  were  between  -1  and  -2  (See  Fig.  3  and  caption  for  details). 
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SWADE  Data  from  Drerinan  et  al..  1996 
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CO  ARE  Data  from  Soloviev  and  Lukas,  2003 


Figure  3.  One  dimensional  turbulence  model  run  results  overlain  with  respective  data  sets  (for  all 
panels  thex-axis  is  the  scaled  TKE  dissipation  rate,  eH/F,  and  the  y-axis  is  the  scaled  depth,  z/Hs). 
In  all  cases  the  green  and  red  curves  are  for  runs  with  injection  of  TKE  flux,  au*w3,  with  a  equal  to 
150  and  100  respectively.  All  TKE  injections  are  at  a  depth  of  Hs/2  except  for  the  purple  curve  on 
panel  g  where  we  used  Hs/4.  Observational  data  were  scaled  with  F,  when  available,  or  by  cai*J, 
with  a  equal  to  150  and  100,  except  for  the  COARE  data  (panel  g)  where  we  used  a  =  50;  the  black 
heavy  line  represents  the  wave  scaling  as  suggested  by  Terray  et  al,  1996.  Note:  model  runs  include 
both  the  shear-free  wave  effects  and  the  shear  generated  TKE  via  the  “ regular  ”  wall-layer 

mechanism. 


RESULTS 

Firstly,  our  results  indicate  that  for  the  majority  of  the  data  sets  analyzed  here  the  widely  used  wall 
layer  parameterization  severely  underestimates  TKE  dissipation  rates  in  the  layer  beneath  wind -waves. 

•5 

In  most  cases  we  found  e/(u*  /kz)  »  1  down  to  depths  of  several  wave  heights,  but  usually  confined  to 

less  than  10  m  (exceptions  are  the  observations  of  Gregg,  1987,  when  s/(u*  /kz)  »  1  down  to  depths 

•2 

of  25  m).  Average  dissipation  estimates  regularly  exceeded  u*  /kz  by  at  least  a  factor  of  4  while  most 

■j 

where  larger  by  an  order  of  magnitude  or  so  (largest  value  was  more  30  times  u*  /kz). 

The  wave  parameterization  proposed  by  Terray  et  al.  (1996)  was  found  to  hold  for  several  of  the  data 
sets  in  which  enhanced  dissipation  rates  were  observed  (Fig.  3).  Specific  examples  are  the  data  sets 
collected  from  a  submarine  (Osborn  et  al.,  1992),  a  free  rising  profiler  (Anis  and  Mourn,  1995),  a 
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quasi-horizontal  glider  (Greenan  et  ah,  2001),  as  well  as  part  of  the  data  collected  with  the  free  falling 
Epsonde  profiler  (Greenan  et  ah,  2001).  An  exception  is  data  collected  with  the  glider  and  Epsonde 
when  swell  was  observed  in  addition  to  wind-waves.  Another  exception  are  the  data  from  COARE 
(Fig.2,  upper  left  panel;  Soloviev  and  Lukas,  2003);  Although  no  information  regarding  swell  was 
available  for  COARE  data,  these  data  seem  to  mimic  the  behavior  of  data  from  Greenan  et  al.  (2001) 
when  swell  was  present,  and  lie  consistently  above  the  best  fit  of  the  WAVES  data:  sHs/F  =  0.3(z/Hsf 
Another  reason  for  this  apparent  discrepancy  might  be  that  a  result  of  underestimating  F,  the  rate  of 
energy  input  from  the  wind  to  the  waves.  Here  we  have  estimated  F  from  the  surface  friction  velocity 
and  the  significant  wave  height,  which  were  the  only  available  parameters.  But,  when  we  used  F  =  50 
u*J  to  scale  the  observed  COARE  dissipation  rates,  the  data  seems  to  collapse  around  the  scaling 
proposed  by  Terray  et  al.  (1996)  (Fig.  3g). 

The  one  dimensional  model  runs  using  the  GOTM  turbulence  model  with  the  generic  two-equation 
model  seem  to  provide  reasonable  estimates  of  both  the  vertical  rate  of  decay  of  dissipation  as  well  as 
the  right  orders  of  magnitude.  A  transition  to  wall-layer  behavior  occurs  at  a  depth  of  ~10 Hs.  We  note 
that  the  runs  that  closely  resemble  the  observed  data  exhibit  a  vertical  decay  rate  of  ~z~2'5  and  a  TKE 
decay  rate  of  ~ z '  (not  shown).  Injection  depths  of  TKE  flux  of  z/Hs  =  0.5  were  found  to  work  well  in 
most  cases,  except  for  the  COARE  data  where  z/Hs  =  0.25  was  found  to  produce  better  results.  The 
sensitivity  of  various  models  and  scaling  to  this  parameter  (in  some  cases  referred  to  as  the  roughness 
length,  z0)  has  already  been  discussed  in  the  literature. 

Estimates  of  the  fraction  of  EI0  dissipated  in  the  OBL  suggested  by  Oakey  and  Elliott  (1982),  i.e.  ~1%, 
was  found  to  be  an  underestimate  when  SL  dissipation  values  were  taken  into  consideration.  Most  of 
the  data  sets  examined  show  that  5%  and  up  to  more  than  10%  of  E10  are  in  fact  dissipated  in  the  OBL. 
This  conforms  closely  to  the  values  suggested  by  Richman  and  Garrett  (1977).  Some  of  the  data  we 
have  analyzed  suggest  that  when  the  waves  are  fetch  limited  and/or  relatively  young,  there  might  be  a 
better  conformity  to  features  predicted  for  the  wall  layer  (e.g  the  fraction  of  Ew  dissipated  in  the  OBL 
is  closer  to  the  predicted  1%).  This  behavior  is  consistent  with  other  evidence  (e.g.  Thorpe,  1992)  and 
we  are  currently  examining  the  effect  of  wave  age  more  closely. 

IMPACT/APPLICATION 

Results  of  this  work  will  improve  TKE  dissipation  parameterization  schemes  used  in  oceanic  models 
and  our  understanding  of  turbulence  and  small-scale  processes  in  the  oceanic  near  surface  layer  in  the 
presence  of  waves. 

TRANSITIONS 

Results  are  currently  discussed  with  several  other  investigators  in  our  research  group. 

A  poster  on  “Statistics  of  Overturns,  APEF,  and  TKE  Dissipation  in  Stratified  Flows”  has  been  given 
at  the  AGU  fall  meeting,  Dec.  2002. 

A  presentation  on  “TKE  Dissipation  Rates  in  the  Aquatic  Wave  Zone”  will  be  given  at  the  Ocean 
Sciences  Meeting  in  Jan.  2003. 

RELATED  PROJECTS 

Other  CBLAST  projects. 


7 


REFERENCES 


Anis,  A.,  and  J.  N.  Mourn,  1995:  Surface  wave-turbulence  interactions:  scaling  s (z)  near  the  sea 
surface.  J.  Phys.  Oceanogr.,  25,  2025-2045. 

Crawford,  G.  B.,  and  W.  G.  Large,  1996:  A  numerical  investigation  of  resonant  inertial  response  of  the 
ocean  to  wind  forcing.  J.  Phys.  Oceanogr.,  26,  873-891. 

Drennan,  W.  M.,  M.  A.  Donelan,  E.  A.  Terray,  and  K.  B.  Katsaros,  1996:  Oceanic  turbulence 
dissipation  measurements  in  SWADE.  J.  Phys.  Oceanogr.,  26,  808-815. 

Greenan,  B.  J.  W.,  N.  S.  Oakey,  and  F.  W.  Dobson,  2001:  Estimates  of  dissipation  in  the  ocean  mixed 
layer  using  a  quasi  horizontal  microstructure  profder.  J.  Phys.  Oceanogr.,  31,  992-1004. 

Gregg,  M.  C.,  1987:  Structures  and  fluxes  in  a  deep  convecting  mixed  layer.  Dynamics  of  the 

Ocean  Mixed  Layer,  P.  Muller  and  D.  Henderson,  Eds.,  Hawaii  Institute  of  Geophysics  Special 
Publication,  1-23. 

Kitaigorodskii,  S.  A.,  M.  A.  Donelan,  J.  L.  Lumley,  and  E.  A.  Terray,  1983:  Wave-turbulence 
interaction  in  the  upper  ocean.  Part  2:  Statistical  characteristics  of  wave  and  turbulent 
components  of  the  random  velocity  field  in  the  marine  surface  layer.  J.  Phvs.  Oceanogr.,  13, 
1988-1999. 

Oakey,  N.  S.,  and  J.  A.  Elliott,  1982:  Dissipation  within  the  surface  mixed  layer.  J.  Phvs.  Oceanogr., 

12,  171-185. 

Osborn,  T.,  D.  M.  Farmer,  S.  Vagle,  S.  A.  Thorpe,  and  M.  Cure,  1992:  Measurements  of  bubble 
plumes  and  turbulence  from  a  submarine.  Atmos. -Ocean,  30,  419-440. 

Richman,  J.,  and  C.  Garrett,  1977:  The  Transfer  of  Energy  and  Momentum  by  the  Wind  to  the  Surface 
Mixed  Layer.  J.  Phys.  Oceanogr.,  7,  876-881. 

Soloviev,  A.,  and  R.  Lukas,  2003:  Observation  of  Wave  Enhanced  Turbulence  in  the  Near 

Surface  Layer  of  the  Ocean  during  TOGA  COARE.  Deep-Sea  Research  I,  50,  371-395. 

Stewart,  R.  W.,  and  H.  L.  Grant,  1962:  Determination  of  the  rate  of  dissipation  of  turbulent  energy  near 
the  sea  surface  in  the  presence  of  waves.  J.  Geophys.  Res.,  67,  3177-3180. 

Terray,  E.  A.,  M.  A.  Donelan,  Y.  C.  Agrawal,  W.  M.  Drennan,  K.  K  Kahma,  A.  J.  Williams  III,  P.  A. 
Hwang,  and  S.  A.  Kitaigorodskii,  1996:  Estimates  of  kinetic  energy  dissipation  under  breaking 
waves.  J.  Phys.  Oceanogr.,  26,  792-807. 

Terray,  E.  A.,  M.  A.  Donelan,  Y.  C.  Agrawal,  W.  M.  Drennan,  K.  K  Kahma,  A.  J.  Williams  III,  P.  A. 
Hwang,  and  S.  A.  Kitaigorodskii,  1997:  Reply.  J.  Phys.  Oceanogr.,  27,  2308-2309. 

Thorpe,  S.  A.,  1992:  Bubble  clouds  and  the  dynamics  of  the  upper  ocean.  Q.  J.  R.  Meteorol. 
Soc.,  118,1-22. 

Umlauf,  L.,  and  H.  Burchard,  2003:  A  generic  length-scale  equation  for  geophysical  turbulence.  J. 

Mar.  Res.,  61,  235-265. 


8 


